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Abstract. We present a new method to fit the variations of both coordinates of a VLBI component as a function 
of time, assuming that the nucleus of the radio source contains a binary black hole system (BBH system). The 
presence of a BBH system produces 2 perturbations of the trajectory of the ejected VLBI components. By using 
only the VLBI coordinates, the problem we have to solve reduces to an astrometric problem. Knowledge of the 
variations of the VLBI coordinates as a function of time contains the kinematical information, thus we are able 
to deduce the inclination angle of the source and the bulk Lorentz factor of the ejected component. Generally, 
there is a family of the BBH system producing the same fit to our data. To illustrate this method, we apply it 
to the source 1807+784. We find that the inclination of the source is io = 5.8° g and the VLBI component is 
ejected with a bulk Lorentz factor of 7 = 3.7 Ig'j ■ We determine the family of the BBH system which provides 
the best fit, assuming at first that the masses of the 2 black holes are equal and then that the masses are different. 
Each family of BBH systems is characterized by Tp/Tt ~ 1.967, where Tp and Tt are the precession period of the 
accretion disk of the black hole ejecting the VLBI component and the orbiting period of the BBH system. 

Key words. Astrometry; BL Lacertae objects: individual: 1803+784; Galaxies: jets 



• ^ 1. Introduction 1. The precession of the accretion disk, 

' riTTTTn 2- motion of the two black holes around the gravity 

H ; previous articles, (Britzen et al. [M and Lobanov & ^^^^^^^ ggjj ^y^^^^ 

5t 1 Roland l2005p we have shown that VLBI and optical obser- 
vations of compact radio sources can be explained if their In addition to these two perturbations, the ejected flow 
nuclei contain a binary black hole system (BBH system), is perturbed by the slow motion of the BBH system around 
However, for most of the compact radio sources we have the center of gravity of the galaxy. This third perturba- 
only VLBI observations. tion is often observed in compact radio sources and is re- 
in this article we present a method to fit the VLBI sponsible for the slow bending of the VLBI jet, i.e., the 
observations using a BBH system. mean ejection direction changes slowly as the distance of 
This method requires knowledge of the variations of ^^'^ component from the core increases. This bending be- 
the two coordinates of the VLBI component as a function comes prominent at a distance from the core of a few mil- 
of time. As these observations contain the kinematical in- Harcsecond (mas). Thus, by studying the motion of the 
formation needed, we are able to deduce the inchnation VLBI component in the innermost part of the VLBI jet, 
angle of the source and the bulk Lorentz factor of the can expect that the influence of the slow motion of 
ejected component BBH system is negligible compared to the perturba- 
We propose a geometric model that assumes that the ^ioi^^ induced by the precession of the accretion disk and 
motion of a component ejected by the nucleus of the radio ^he motion of the black holes around the gravity center of 
source is perturbed by two different motions, namely: BBH system. The inclusion of this slow motion m this 

model would be difficult as the characteristics of the BBH 
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of a component moving in the first mas of the VLBf jet 
are needed. Using just the VLBI coordinates and a geo- 
metric model, the problem we have to solve reduces to an 
astrometric one. 

The method presented in Lobanov & Roland (2005) 
was not consistent. In this article we built a consistent 
method that solves the problems found in Lobanov & 
Roland (2005). 

In Section 2 we will present the details of the model 
and in Section 3 the new method. 

To illustrate the method, we apply it to the source 
S5 1803+784 and we show how, from the knowledge of 
the variations of the coordinates X{t) and Y{t) of a VLBI 
component, we can derive 

— the inclination angle of the source, 

— the bulk Lorentz factor of the ejected component, 

— the angle between the accretion disk and the plane of 
the BBH system, 

— the size of the BBH system, 

— the ratio Tp/T),, where Tp is the precession period of 
the accretion disk and the orbiting period of the 
BBH system, 

— the origin of the ejection of the VLBI component, 

— the duration of the ejection of the plasma responsible 
for the VLBI component. 

2. The model 

2.1. Introduction: The two-fluid model 

We will describe the ejection of a VLBI component in the 
framework of the two-fluid model (Sol et al. 1989, Pelletier 
& Roland 1989, 1990, Pelletier & Sol 1992). The two-fluid 
description of the outflow is adopted with the following 
assumptions : 

1. The outflow consists of an e~ —p plasma (hereafter the 
jet) moving at mildly relativistic speed Vj < 0.4 x c and 
an plasma (hereafter the beam) moving at highly 
relativistic speed (with corresponding Lorentz factor 
7fc < 30). 

2. The magnetic field lines are parallel to the flow in the 
beam and the mixing layer, and are toroidal in the jet 
(see Figure 1 of Lobanov & Roland 2005). 

The e~ — p jet carries most of the mass and the kinetic 
energy ejected by the nucleus. It is responsible for the for- 
mation of kpc-jets, hot spots and extended lobes (Muxlow 
et al. 1988, Roland et al. 1988, Roland & Hetem 1996). 
The relativistic beam moves in a channel through the 
mildly relativistic jet and is responsible for the formation 
of superluminal sources and their 7-ray emission (Roland 
et al. 1994). The relativistic beam can propagate if the 
magnetic field B is parallel to the flow in the beam and 
in the mixing layer between the beam and the jet and 
if it is greater than a critical value (Pelletier et al. 1988, 
Achatz & Schlickeiser 1993). The magnetic field in the jet 
becomes rapidly toroidal (Pelletier & Roland 1990). 



The observational evidence for the two-fluid model has 
been discussed by e.g. Roland & Hetem (1996). Recent 
observational evidence for relativistic ejection of an 
beam come from the 7-ray observations of MeV sources 
(Roland & Hermsen 1995, Skibo et al. 1997) and from 
VLBI polarization observations (Attridge et al. 1999). 

The formation of X-ray and 7-ray spectra, assuming 
relativistic ejection of beams, has been investigated by 
Marcowith et al. (1995, 1998) in the case of Centaurus A. 

The possible existence of VLBI components with two 
different speeds has been recently pointed out in the case 
of the radio galaxies Centaurus A (Tingay et al. 1998), 
Virgo A (Biretta et al. 1999) and 3C 120 (Gomez et 
al. 2001). If the relativistic beam transfers some energy 
and/or relativistic particles to the jet, the relativistic par- 
ticles in the jet will radiate and a new VLBI component 
with a mildly relativistic speed will be observed (3C 120 
is a good example of a source showing this effect). 

2.2. The geometry of the model 

We will call the angle between the accretion disk and the 
orbital plane {XOY) of the BBH system. The ejection of 
the component will be on a cone with its axis in the Z'OZ 
direction and opening angle fl. We will assume that the 
line of sight is in the plane (YOZ) and forms an angle io 
with the axis Z'OZ (see Figure [1]). The plane perpendic- 
ular to the line of sight is the plane (r/OX). We will call 
A2 the angle of the rotation in the plane perpendicular to 
the line of sight, in order to transform the coordinates r/ 
and X into coordinates N (North) and W (West), which 
are directly comparable with the VLBI observations. We 
have 

W = X cos(AS) — [z sin(io) + y cos(io)) sin(AS) , (1) 



N = X sin(A^) -|- {z sin(io) -t- y cos(io)) cos(AE!) . (2) 



2.3. The precession model 

Here we describe the precession of the accretion disk. 

The coordinates of a component moving in the per- 
turbed beam are given by 

Xc = Ro{z)cos{uOpt - kpz{t) + (po) , (3) 



j/c = Ro{z)sin{ujpt - kpz{t) + (f>o) , (4) 



Zc = Zc{t) , (5) 
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Fig. 1. The geometry of the problem 

where uip = 2-K/Tp, Tp is the precession period, and kp is 
defined by 



kp = 2Tr /TpVa , 



(6) 



where T4 is the speed of the propagation of the perturba- 
tions and a free parameter of the problem. 

Wc will assume that the amplitude of the perturba- 
tion first increases linearly, and we take the form of the 
amplitude R{zc{t)) to be 



R{Zc{t)) 



RgZeit) 

{a + zS)) ' 
where a is 

a = Ro/{2 tanfl) . 
2.4. The binary system model 



(7) 



(8) 



To explain the origin of the precession of the accretion 
disk, wc propose that the nucleus hosts a BBH system. 

As stated above, the two black holes orbit in the plane 
(XOY), and the origin of our coordinate system is centered 
on the mass center of the system. The elliptical orbit is 
given by 



r = 



P 



1 



(9) 



where e and p are respectively the eccentricity and the 
parameter or the semi-latus rectum of the orbit. 

We will assume that the two black holes have circular 
orbits, i.e. e = and wo will define the black hole that 
ejects the VLBI component with index 1. Its coordinates 
are : 



Yi{t) 



Mo 



M1+M2 



psm{ip{t)) . 



(11) 



As the orbits are circular, we have ijj{t) = Uht + 'Re- 
writing 



xi = yi 
we have 



pM2 



Ml +M2 



Mo 



M1 + M2 



^G(Mi + M2) 



1/3 



(12) 



(13) 



where T(, is the period of the BBH system. 

We define R^bh the distance between the two black 
holes as the size of the BBH system, it is: 



R 



■bbh 



^G{Mi + M2) 



1/3 



(14) 



The coordinates of black hole 1 can be written 



Xi{t) = xi cos{uibt + V'o) , 



Yi{t) = yi sm{u)bt + V'o) • 



(15) 



(16) 



Finally, we suppose that the perturbation of the beam 
is damped with a characteristic time Td- 

For VLBI observations the origin of the coordinates is 
black hole 1, i.e. the black hole ejecting the VLBI com- 
ponents. Therefore, the coordinates of the moving com- 
ponents in the frame of reference where black hole 1 is 
considered the origin are 



Xc = [Ro{z) cos{ajpt - kpz{t) + 4>o) 

+ XiCOs{u)bt — kbz{t) + ipo) — X\ cos(V'o)] 
exp{-t/Td) , (17) 



yc = [Ro{z)s\n{i^pt - kpZ{t) + (t)o) 

+ yi smiujbt - hz{t) + ipo) - yi sin(Ro)] 
exp{-t/Ta) , (18) 



where coh = 27r/T(„ and kh is defined by 



(19) 



X,{t) = - 



Mo 



Ml +M2 



pcos{tp{t)) 



(10) 



kb = 



27r 

TbVa' 



(20) 
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The differential equation governing the evolution of 
Zc{t) can be obtained through the relation for the speed 
of the component 



dt 



dt 



\ dt 



(21) 



where is related to the bulk Lorentz factor by Vc/c — 
V(l-l/7?). 

Using (O, (g]) and dS]), we find from ^ that dz^dt 
is the solution of the equation 



A 



dzc 
dt 



. 



(22) 



A calculation of the coefficients A, B and C can be 
found in Appendix I. 

Equation admits two solutions corresponding to 
the jet and the counter-jet. 

We assumed that the line of sight is in the plane (YOZ) 
and makes an angle io with the z axis (see Figured]). Thus 
following Camenzind & Krockenberger (1992), if we call 
the angle between the velocity of the component and the 
line of sight we have 



cos{e{t)) 



dye . . dzc . , , 

— lo + "^cos lo ] /Vc 



(23) 



The Doppler beaming factor S, characterizing the 
anisotropic emission of the moving component, is 



Sc{t) 



1 



7c [1 - PccMm)] ' 

where Pc = Vc/c. The observed flux density is 



(24) 



5c = ^<5c(t)'+"'-(l + z)^-"'- 



JcdV , 



(25) 



where D is the luminosity distance of the source, z its 
redshift, jc is the emissivity of the component, and is 
the synchrotron spectral index (a negative definition of 
the spectral index, S oc v^" is used). As the component is 
moving relativistically toward the observer, the observed 
time is shortened and is given by 



tobs 



[1- f]cCOs{0{t' j)]{l+z)dt' 



(26) 



3. The global method 

3.1. Introduction 

In Lobanov & Roland (2005) we provide a method to de- 
termine the characteristic parameters of the BBH system 
using radio and optical observations. The method consists 
of two steps, i.e. in a first step we use VLBI observations 
to model the precession (without a BBH system) and in 



a second step we use optical observations to obtain the 
characteristic parameters of the BBH system. 

The above described method has the following prob- 
lems: 

1. Using a simple precession model, we find for 3C 345 
that a bulk Lorentz factor increasing with time is nec- 
essary. 

2. Using the parameters of the precession found in the 
first step, the BBH solution obtained in the second 
step is not necessarily consistent with the precession 
solution found previously. Indeed, in the limit M2 — > 
the BBH solution is not necessarily able to reproduce 
the results of the precession model found in the first 
step. The limit M2 corresponds to a single black 
hole and the precession of the accretion disk is due to 
the Lens-Thirring effect in that case. 

The above problems can be solved if we directly model 
the VLBI observations with a BBH system instead of a 
simple precession model. In the BBH system model, the 
bulk Lorentz factor is constant, and the model explains 
the apparent variations of the speed of the VLBI compo- 
nent when it escapes from the nucleus (see Figure [SJ the 
apparent speed of the ejected component changes by a fac- 
tor of four with a constant bulk Lorentz factor and it is 
not necessary to involve any acceleration or decceleration 
of the VLBI component). Since the BBH system and the 
precession solution are obtained simultaneously, they are 
obviously self-consistent. 

We call the method we present in this article the global 
method. 

In the global method, we calculate the projected tra- 
jectory on the plane of the sky of a component ejected 
by a BBH system and we determine the parameters of 
the model to simultaneously produce the best fit for both 
the West and North coordinates, i.e. Wc{t), Nc{t). So, the 
parameters found are such that minimize 

x't=x'{WS))+x'{N,{t)) , (27) 

where x^(Wc(i)) and x^(-^c(0) are the calculated by 
comparing the VLBI observations with the calculated co- 
ordinates Wc{t) and Nc{t) of the component. 

3.2. The coordinates of the VLBI component 



Solving (im, we determine the coordinate Zc{t) of a point 
source component ejected relativistically in the perturbed 
beam. Then, using (|17p and (fT8|) . we can find the coordi- 
nates Xc{t) and yc{t) of the component. In addition, for 
each point of the trajectory, we can calculate the deriva- 
tives dxc/dt, dyc/dt, dzc/dt and then deduce cos 6 from 
(pa . Sc from dm), from ^ and tobs from (^5)1 . 

When the coordinates Xc{t), yc{t) and Zc{t) have been 
calculated, they can be transformed to wdt) (West) and 
nc{t) (North) coordinates using ^ and ([2|). 

As explained in Britzen et al. (|2001|) and Lobanov & 
Roland (2005), the radio VLBI component has to be de- 
scribed as an extended component along the beam. Let us 
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call rirad the number of points (or steps along the beam) 
for which we integrate, in order to model the component. 
The coordinates Wc(i), Nc{t) of the VLBI component are 
then 

Wc{t) = {^^W,,{t)^ /Urad , (28) 

Nc{t) = Z"™'^ • (29) 

and can be compared with the observed coordinates of the 
VLBI component. 

3.3. The parameters of the model 

In this section, we list what a priori the free parameters 
of the model are: 

— io the inclination angle, 

— (po the phase of the precession at i = 0, 

— AS the rotation angle in the plane perpendicular to 
the line of sight (see H]) and 

— Q. the opening angle of the precession cone (see ([8|)), 

— Ro the maximum amplitude of the perturbation (see 

. . 

— Tp the precession period of the accretion disk, 

— Tfi the characteristic time for the damping of the beam 
perturbation, 

— Ml the mass of the black hole ejecting the radio jet, 

— M2 the mass of the secondary black hole, 

— 7c the bulk Lorentz factor of the VLBI component, 

— V'o the phase of the BBH system at t = 0, 

— Tb the period of the BBH system, 

— to the origin of the ejection of the VLBI component, 

— Va the propagation speed of the perturbations, 

^ TT-rad IS the number of steps to describe the extension 
of the VLBI component along the beam. 

To begin with, we assume that Mi = M2 and when 
the corresponding solution is obtained, we calculate the 
family Mi(M2) which provides the same fit. 

So, the problem we have to solve is a 14 free parameters 
problem. 

If, in addition to the radio, optical observations are 
available that peak in the light curve, this optical emission 
can be modelled as the synchrotron emission of a point 
source ejected in the perturbed beam (Britzen & al. 120011 
and Lobanov & Roland 2005). This short burst of very 
energetic relativistic is followed immediately by a very 
long burst of less energetic relativistic . This long burst 
is modelled as an extended structure along the beam and 
is responsible for the VLBI radio emission. In that case 
the origin to of the VLBI component is the beginning of 
the first peak of the optical light curve and is not a free 
parameter of the model. 

We have to investigate the different possible scenar- 
ios with regard to the sense of the rotation of the accre- 
tion disk and the sense of the orbital rotation of the BBH 



system. These possibilities correspond to ± ujp{t — z/Va) 
and ± ijjf,(i — z/Va). As the sense of the precession is al- 
ways opposite to the sense of the orbital motion, we will 
study the two cases denoted by H — and — (- where we have 

UJp{t-z/Va), -LObit-z/Va) and -UJp{t-z/Va), UJb{t-z/Va) 

respectively. 

3.4. The method to solve the problem 

To find a solution for the above described problem we use 
the following method. 

As mentioned before we start with the assumption that 
Ml — M2, i.e. that the two masses of the BBH system are 
equal. 

First, we find the inclination angle that provides the 
best fit. To do that we minimize Xti^o) (see (P7))) when 
the inclination angle varies gradually between two values. 
At each step of io, we determine each free parameter A 
such that dxt /dX = and Xt presents a minimum. 

Furthermore, using the inclination angle determined 
previously, we explore the space of the solutions, for a 
varying mass of the BBH system (while still assuming 
Ml = M2). This allow us to find whether the solution 
of the BBH system presents a degeneration or if there are 
other solutions, with different masses, that fit the obser- 
vations. When exploring the solutions space, we always 
vary one parameter in a step-wise manner, with each step 
minimizing for each of the free parameters. 

The space of the solutions can be explored for each of 
the free parameters if necessary. 

Because the problem is a non-linear one, we calculate 
again the variations of Xt (*o) for the best solution found 
previously, starting from the inclination angle obtained in 
the first step. Where Xt (*o) reaches its minimum, we have 

I o • I ^ ^ (*o ^ io,mvn) i (30) 

V / mm 

and the 1 a error bar, (Aio)io^, corresponding to the pa- 
rameter io is then given by 

{^io)ia = l/A . (31) 

This assumes that around the minimum, Xt(*o) is a 
parabola, however, for large variations of io, the parabola 
approximation is not valid and a better determination of 
the 1 a error bar can be obtained using the definition 

{Aio)la = |jo(xL« + 1) - *o(xLrJI , (32) 

which provides two values {Aio)ia+ and (Aio)i(j- (see 
Lampton et al. 1976 and Hebrard & al. 2002). 

Because we calculated Xt (*o) by minimizing XtW 
each step and for each free parameter A, we can deduce the 
range of values corresponding to 1 cr for each parameter 
when the inclination varies between io = io,min — {A.io)ia-- 
and io = io,min + (Aio)icr+ • 
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Finally, when the best solution corresponding to AIi ~ 
M2 is obtained, we can determine the family of BBH sys- 
tems with Ml ^ M2 that provides the same fit. 

One important point is that using this phenomeno- 
logical method, we do not have the proof that the mini- 
mum found is unique, since for a completely different set 
of values for the parameters of the problem, another mini- 
mum could exist. However, as we will explore a wide range 
of inclination angles, i.e. 1° < io < 25° and a equally 
wide range of BBH system masses, i.e. 10^ x Mq < M < 
10^^ X Mq, we minimize the possibility of missing the 
best solution. Another way to overcome this difficulty is 
to explore the space of possible values of the parameters 
using, for instance, a Monte Carlo Markov chain algorithm 
(MCMC algorithm). However this is out of the scope of 
this article. 

4. Application to S5 1803+784 

4.1. The radio morphology of 55 1803+ 784 

The blazar S5 1803-^784 (z « 0.68, Lawrence et al. 1987, 
Stickel et al. 1993) is an intraday-variable (IDV) source 
with rapid flux-density variations in the optical and radio 
regime (Wagner & Witzel 1995) on timescales as short 
as 50 minutes in the optical (Wagner et al. 1990). The 
source has been observed and studied with different VLBI 
arrays at a wide range of frequencies (e.g., Eckart et al. 
1986, 1987; Witzel et al. 1988; Chariot 1990; Strom & 
Biermann 1991; Fey et al. 1996; Gabuzda 1999; Gabuzda 
& Cawthorne 2000; Gabuzda & Chernetskii 2003; Ros et 
al. 2000, 2001; Britzen et al. 2005a, 2005b). Global VLBI 
observations reveal the pc-scale jet of S5 1803-f 784 to be 
oriented in the East- West direction (e.g., Britzen et al. 
2005b), while the large-scale structure comprises a dom- 
inant core component and a weak secondary component 
~ 45" away from the core at position angle « —166° 
south-south-west (Antonucci et al. 1986; Britzen et al. 
2005a) and makes this source a "misaligned" object (e.g., 
Pearson & Readhead 1988). Strom & Biermann (1991) 
also detected this secondary component, as well as a 
bridge connecting it with the core, at 1.5 GHz obser- 
vations with the Westerbork Synthesis Radio Telescope 
(WSRT). World-array VLBI observations revealed that 
the jet bends by around 90° at a core separation of about 
0.5 arcseconds towards the South (Britzen et al. 2005a). 
The Westerbork Synthesis Radio Telescope (WSRT) ob- 
servations resolved the bridge connecting the core with 
the secondary component some 45" to the south. Wiggles 
in the ridge of this emission suggest that, as on the 100- 
pc scale, the jet may also oscillate on 50-100 kpc scales. 
In addition, amorphous emission to the north, as well as 
an extended, halo-like component around the nucleus has 
been detected (Britzen et al. 2005a). 

On pc-scales the source shows a pronounced jet, with 
prominent jet components located at distances between 
1.4 and 12 mas from the core (e.g., Eckart et al. 1986). 
Geodetic and astronomical VLBI data obtained at var- 



ious epochs between 1979 and 1985 indicated that the 
component at 1.4 mas is stationary (e.g., Schalinski et 
al. 1988; Witzel et al. 1988). Britzen et al. (2005b) dis- 
cuss significant position shifts for the 1.4 mas component 
with displacements between r ~ 0.7 mas and r 1.5 
mas. Obviously the improved time sampling of the geode- 
tic VLBI data led to the detection of systematic position 
variations for this component regarded as stationary on 
the basis of the less frequently targeted astronomical ob- 
servations. This oscillatory behaviour is explained as the 
result of a reconfinement shock (Britzen et al. 2005b). A 
recent analysis of 94 epochs of VLBI observations obtained 
at V = 1.6, 2.3, 5, 8 15, and 43 GHz reveals that that the 
jet structure within 12 mas from the core can most easily 
be described by four jet components that remain at similar 
core separations in addition to the already known "oscil- 
lating" jet feature at ~ 1.4 mas. We show the pc-scale 
structure of 1803-1-784 in an image obtained in VLBA ob- 
servations performed at 15 GHz (April 2004) in Fig. [2l 
The jet components are indicated. In addition to these 
"stationary" components, we find one much fainter com- 
ponent, B3, moving with apparent superluminal velocities. 
This component most likely coincides with the component 
seen in 43 GHz maps by Jorstad et al. (2005). However, at 
15 GHz this component is much fainter compared to the 
other components. 



Clean I map. Array: BFHKLMNOPS 
1805 + ^84 at 15-365 GHz 2004 Apr 10 




-5 -10 -1 

Right Ascension (mas) 
Map center: RA: 18 00 45-084, Dec: +78 28 04.019 [2000.0) 
Mop peak: 1.26 Jy/beam 

Contours %: -0.12 0.12 0-24 0.48 0.96 1.92 3.84 

Contours %: 7.68 15.4 30.7 61.4 

Boom FWHM: 0-6 x 0.579 (mos) at -84.5° 



Fig. 2. The pc-scale structure of 1803-1-784 as obtained in 
15 GHz VLBA observations in April 2004. The jet com- 
ponents up to a core separation of ^12 mas are marked. 

Evidence is presented for quasi-periodicities in the 
variability of the core separation, the position angle and 
the flux of the "oscillating" jet components on timescales 
comparable to those derived from the total flux density 
changes (Britzen et al., in prep.; Kudryavtseva et al., in 
prep.). Superluminal motion has been detected in the jet 
of S5 1803-f 784 at 43, 22, (Krichbaum et al. 1993) and 
8.4 GHz (Britzen et al. 2005b; Britzen et al., in prep.). 
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In Britzen et al. (2005b) we find and discuss that three 
jet components approach the brightest and so-called "sta- 
tionary" component (at ~ 1.4 mas at 8.4 GHz) with an 
apparent super luminal motion of 8 - 11 c. In this paper we 
show that the "stationary" component oscillates, - under 
the assumption that this is the brightest jet component in 
each epoch - and we discuss several possible explanations. 
Based on the most recently performed investigations of a 
much larger database (94 epochs of VLBI observations) 
covering a longer time span (of almost 20 years) and sev- 
eral frequencies (1.6, 2.3, 5, 8.4, and 15 GHz) Britzen et 
al. find that the brightest jet components near the core do 
not reveal fast radial motion. Instead, components tend 
to remain at similar core separations but show significant 
motion with regard to the position angle. Some indication 
of apparent superluminal motion between 5 c and 7 c 
for the 25 mas jet component has been derived from 6 and 
18 cm observations (Britzen et al. 2005a). VLBI observa- 
tions at 43 GHz revealed, for the first time, evidence for 
a jet structure described as "helical" by Krichbaum (e.g., 
1990). A curved-jet morphology is found at all scales in- 
vestigated so far (e.g., Britzen et al. 2005a). 

4.2. The observations 

For the present analysis we use the archival VLBI obser- 
vations at two frequencies, 8 and 15 GHz, spanning almost 
20 years of observations (for details see Britzen et al., in 
prep.). S5 1803-1-784 has been observed with VLBI at A 
= 15 GHz by Perez-Torres et al. 2000, Kellermann et al. 
(1998, 2004), and Zensus et al. (2002) between 1994.67 
and 2005.68 in 13 epochs. The 8 GHz observations were 
performed by Ros et al. (2000, 2001) and Perez- Torres et 
al. (2000) from 1986.21 to 1993.95 in 41 epochs. The data 
had been fringe-fitted and calibrated before by the indi- 
vidual observers (for details, see the references given in 
Britzen et al., in prep.). We performed modelfitting of cir- 
cular Gaussian components within the Difmap package. In 
order to find the optimum set of components and param- 
eters, we fitted each data set starting from a point-like 
model. Circular components have been chosen in order 
to simplify the comparison and to avoid unlikely and ex- 
tremely extended elliptical components. For details of the 
data reduction and the modelfit parameters see to Britzen 
et al. (in prep.). 

4.3. Preliminary remarks on the fit 

Before we begin the fit of VLBI components of 
S5 1803-1-784, we indicate some values we will use for the 
cosmological parameters and the error bars. 

The cosmological model we will use is defined by 
Qt = f^A + ^m, with r^A — 0.7, Qm = 0.3 and Ho = 71 
km/s/Mpc for the Hubble constant, fit, flA and flm are 
respectively the total density of the Universe, the density 
of the vacuum and the density of the matter. Calling Zs 



the redshift of the source, the luminosity distance, Di, is 
defined by 

where 



E{zs) = ^0.3(1 + z,f + 0.7(1 + Zs) . (34) 

With Zs « 0.68 we have Di k. c x 0.8995/iJo and the 
angular distance is Da — Di/{1 + z)^. 

The beams corresponding to 8 GHz and 15 GHz ob- 
servations are respectively w 1.0 mas and w 0.5 mas. The 
minimum error bars we will use are Aw = 0.08 mas and 
An — 0.08 mas for the West and North coordinates at 
8 GHz and Aw = 0.04 mas and An = 0.04 mas for the 
West and North coordinates at 15 GHz. 

4.4. VLBI components CO, CI and the trajectory of 
the ejected plasma component 

The jet of S5 1803-f784 consists of five "oscillating" jet 
components, CO at ~ 0.3 mas, CI at ~ 0.7 mas, Ca at 
~ 1.4 mas, C2 at ~ 2 mas and C4 at about ~ 4 mas (for 
the details see Britzen et al. in prep.). In this paper we 
will discuss only the motion of the two innermost com- 
ponents CO and CI. Previous analysis has shown that the 
brightest component Ca, which is quasi-stationary and has 
almost the same core separation, is changing its position 
very slowly with a speed of about 0.04 mas/yr (see Fig. [3]). 
This could mean that all the oscillating jet components of 
S5 1803+784 are moving extremely slowly outwards with 
time. The pertubation propagation speed of the magnetic 
tube (the red curve of Figure [6]) is Va <C c and is the same 
for the magnetic loops corresponding to CO, CI and CA. 
The source S5 1803-1-784 has been observed over a time 
span of 20 years. Thus, we can trace this slow motion of 
the components with time. The distance from the core of 
the component CO is ^ 0.3 mas and given a speed of about 

0. 04 mas/yr it will change its position by ~ 0.7 mas in 10 
years, which is actually the position of the component CI. 
Therefore, we can assume that the component CO, which 
is observed as CO during the first six years of observa- 
tions, (between 1986 and 1992), becomes component CI 
and from epoch 1996.38 follows the trajectory of this com- 
ponent. So, by using the 20 years of VLBI observations of 
S5 1803+784 as follows: 

1. from 1986.21 to 1991.94, VLBI data corresponding to 
CO (8 GHz), 

2. for 1993.95, the VLBI point located between CO and 
CI (8 GHz), and 

3. from 1996.38 to 2005.68, VLBI data corresponding to 
CI (15 GHz), 

we can calculate the trajectory that this new plasma com- 
ponent follows and define the moment of ejection from the 
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core, which we find to be around epoch 1984.5 (see the 
Figures S] , [5] and O . Since there have not been any de- 
tected frequency-dependent effects in the position of the 
jet components (Britzen et al. in prep.), we used 8 and 15 
GHz data together for the derivation of the trajectories of 
the jet components CO and CI. 



0.3 



0.2 
0.1 

1 °° 
> 

-0.1 
-0.2 
-0.3 



1803+784 




• CaBGHz: 1986- 1994 yr5. 
O Ca 1 5 GHz: 1 995 - 2005 yrs. 
A Ca 8 GHz literature data: 1995 - 2001 yrs. 



0.2 0.4 0.6 O.i 



1.0 1.2 1.4 1.6 1.8 2.0 
X [mas] 



Fig. 3. The movement of the Ca component. The red cir- 
cles represent the position in rectangular coordinates of 
the Ca component at 8 GHz during the period 1986 - 
1994. The blue triangles show the position at 8 GHz dur- 
ing 1995 - 2001 and the grey diamonds the position at 15 
GHz from 1995 to 2005. The component moved 0.4 mas 
during ten years of observations. 
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Observed Time 



Fig. 4. The West coordinate of 55 1803 + 784. The small- 
est error bars for 8 GHz observations are 0.08 mas and 
for 15 GHz observations are 0.04 mas 

Around 1993, 9 years after the ejection of the previ- 
ous plasma component, a new plasma component arises 
and follows the magnetic tube in the beam, producing 
a new VLBI component corresponding to CO again (see 
Figure [2]). In fact, the quasi- stationary component CO 
from Britzen et al. (in prep.) corresponds to the quasi- 
stationary perturbation of the magnetic tube in which the 




1995 2000 

Observed Time 



Fig. 5. The North coordinate of 55 1803-^784. The small- 
est error bars for 8 GHz observations are 0.08 mas and for 
15 GHz observations are 0.04 mas 



relativistic plasma propagates. The magnetic tube defor- 
mation propagates with a speed Va c. 

The existence of a quasi-stationary VLBI component, 
which seems to be oscillating, can be explained by the 
fact that firstly the magnetic perturbation of the beam 
(magnetic loop) is mostly perpendicular to the ejection 
direction (see Figure [6]) and secondly the observed time 
between two ejections of the plasma component from the 
nucleus corresponds to the time for a plasma component 
to cross the abovementioned magnetic loop. 

As the ejection direction is west, the magnetic tube 
deformation propagating with the speed Va <^ c will pro- 
duce a slow motion in the west direction. Moreover, VLBI 
observations of all the components show a slow bending 
to the south after the CO and CI components. This slow 
bending can be an effect of the slow motion of the BBH 
system around the gravity center of the galaxy. 

For the VLBI component Ca at ~ 1.4 mas, which is 
situated just after CI (see Figure [2]) slow motion in the 
south-west direction has been detected over the 20 years 
of observations. The fit of the model to the trajectories of 
CO and CI predicts that this component will be observed 
at the present position of Ca and it will move slowly out- 
wards (see Figure El) . This slow motion is due to both the 
slow motion of the magnetic perturbation and of the BBH 
system around the gravity center of the galaxy. 

Since the information on the absolute position of the 
core is lost in intcrfcrometric observations, we are not able 
to determine whether the core is stationary or is itself 
moving. We will investigate the stationarity of the core in 
phase reference observations. 

As indicated in the introduction, the VLBI core will 
move due to: 

— the motion of the black hole ejecting the VLBI com- 
ponent around the gravity center of the BBH system 
and 
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— the motion of the BBH system around the gravity cen- 
ter of the galaxy. 
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Fig. 6. The trajectory of 55 1803 + 784. In addition to 
components CO and CI, we plotted the points correspond- 
ing to the 20 years of observations of the component Ca. 
The solid line shows the fit of a model to the trajectories 
of CO and CI components. The position of the Ca com- 
ponent lies on the fitted trajectory. This shows the slow 
motion of the beam perturbation which is due to the per- 
turbation speed Va ^ c and the slow motion of the BBH 
system around the gravity center of the galaxy. 



— the parameter AS is maintained constant, i.e. ~ 
271°, 

— the parameter Td is kept as a free parameter but with 
a maximal value Td < Td.max = 1000 yrs (0). 

We calculated Xt (*o) when the inclination angle varies 
between 1° and 25°. The result is shown in Figure [7] 




5 10 15 20 25 

io 

Fig. 7. The variations of Xt (*o)- There is only one mini- 
mum in this interval and it is located at io ~ 5.8°. The fit 
has been made using 66 points. As the model uses 12 free 
parameters the solution is well constrained, as shown by 
the variations of the curve x^(*o) 



4.5. Determination of the minimum ofx^iio) 

We first assume Mi = M2, i.e. that the 2 masses of the 
BBH system are equal. 

From the shape of the trajectory, we find that the pre- 
cession is defined by —ujp{t — z/Va), thus the orbital rota- 
tion has to be ujb{t — z/Va). 

As indicated previously, to find the minimum of x^(*o), 
we minimize Xt (^o) as the inclination angle varies gradu- 
ally between two values. At each step of io, we determine 
each free parameter A such that XtW minimal for A. 

To determine the trajectory, we integrate (fT9|) using a 
step At. By calling ac the size of the ejected component, 
the parameter n^ad, describing the length of the VLBI 
component along the beam, can be related to CTc as follows 

(Tc = rirad sin{io) c At/ (DaTnas) . (35) 

As the inclination varies, we change the integration 
step in order to keep Cc = Cst if rirad = Cst. 

As the 2 components CO and CI are within 1 mas and 
the trajectory is complicated, the direction of the ejection, 
i.e. the parameter A5, and the damping time, Td, of the 
beam perturbations cannot be properly constrained by the 
fit of the CO and CI coordinates. To find the solution 



The inclination angle corresponding to the minimum 
is io ~ 5.8° and the determination of the la using (I32|) 
provides {Aio)ia- ~ 1-8 and (Azo)icr+ « 1.7. Thus we 
have 

lo - 5.8 til ■ (36) 

As we calculated the variations of io step by step min- 
imizing x^ for each variable at each step, we can deduce 
the range of values corresponding to la for each parameter 
when the inclination varies between io = io.min — iAio)ia- 

and io = «o,mm + {Aio)la+- 

The bulk Lorentz factor of the VLBI component pro- 
ducing CO and CI is 

7c = 3.7 t'o.l ■ (37) 

The variations of the apparent speed and the variations 
of the distance from the core of the VLBI component are 
shown in Figure [S] During the times when the motion 
of the component is mainly perpendicular to the mean 

^ The value Td^max = 1000 yrs has been chosen to limit the 
expansion of the perturbation after the third VLBI component 
and values Td > Td.max do not significantly change the value 
of x' 
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ejection direction, as the component moves within a loop 
of the perturbated magnetic tube, its distance from the 
core is mostly constant but its apparent speed is high, i.e. 
Vap ~ 3.7 X c. The periods of time corresponding to these 
phases are between 1987 and 1990 for the first one and 
around 2000 for the second one. 




1985 1990 1995 2000 2005 2010 



Observer Time 

Fig. 8. The variations with time of the apparent speed and 
the distance from the core of the ejected component. We 
see that the distance from the core can be mostly constant 
during certain periods of time. During these periods, the 
apparent speed is high. 



We will discuss the results of the BBH system param- 
eters in the next section. However, here we point out an 
important characteristic of the fit of 1807+784. The bet- 
ter fits are obtained when the ratio Tp/Tb has values close 
to 1, 2 or 3. The corresponding values of Xmm given 
in Table 1 



Table 1 : Values of xLin fo^' different ratios Tp/Tb 



Tp/Tb 


r ■ 


« 1 


47.1 


« 2 


34.8 


« 3 


37.6 



From Table 1, we see that the best fit is obtained when 
Tp/Tb ~ 2.0. The results presented in this article corre- 
spond to this case. 

In Table 2, we give the values of the geometric parame- 
ters characterizing the trajectory of the VLBI component 
producing CO and CI. 



Table 2 : Geometric parameters of the fit of CO and CI 
for io ~ 5.8 



ParametGr 


Value 


umt 


ReiTLRrk 


AS 


271 


o 


Fixed 


fl 


1 09 +''-^^ 

l.Vi) _o.24 


o 




Ro 


581 t'o 


pc 




Td 


1000 


yr 




to 


1984.51 +^-^5 


yr 




^rad 


165 tit 


Steps 




{'^ejec^ obs 


0.96 


yr 


Obs Frame 



The parameter Ro is found to be Ro 3> 1, which in- 
dicates that the increase of the perturbation due to the 
precession is in the linear regime, i.e. Ro oc z (see equa- 
tion [7]). In that case, becomes mostly independent of 
Ro- 

The parameter T^ reaches the maximum allowed value. 
The ejection duration of the plasma responsible for the 
VLBI component in the BBH system frame is 

iTejec)bbh = rirad At « 14.47 yr , (38) 

where At is the integration step. Note that it does not 
depend on the observing frequency. 

We give in Table 2 the value of the duration of the 
ejection of the VLBI component in the observer frame. As 
the time compression factor due to the relativistic ejec- 
tion defined by tobs/t (see equation [26l) is « 0.0668, the 
duration of the ejection of the VLBI component in the 
observer frame is {Tejec)obs ~ 0.96 yr. 

The parameter rirad also characterizes the length of the 
VLBI component along the beam. This length is Cc ~ 0.07 
mas. 

The fit has been mqde using 66 points (33 points for 
each coordinate). As the model uses 12 free parameters 
(indeed AS = Cst and Td — Td^max — Cst), the solution 
is well constrained as shown by Figure [7l 

4.6. Determination of the family of BBH systems 

Starting from the solution found in section 14.51 we fix the 
parameters io = io,min, to = 1984.51, rirad = 165 and 
we study the variations of when the mass Mi varies 
between 10^ Mq and 10^^ Mq (with the assumption Mi = 
M2 still holding true). 

This allows us to find whether, for a given mass, unique 
solutions, or families of solutions exit. For the latter case, 
we find that Xt ~ C'st when the parameter Mi varies, i.e. 
the parameter Mi shows a degeneration. 

We find that there are 2 families corresponding to the 
ratio Tp/Tb ~ 2.0. The results corresponding to these fam- 
ilies are given in Table 3. 

Table 3 : Families of solutions with Tp/Tb ~ 2.0 for 



io « 5.8 


Family 


xtiMi) 


Tp/n 


SI 


« 34.8 


1.967 


S2 


« 35.8 


2.01 



The best fit corresponds to the family SI for which we 
have 
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^ « 1.967 



(39) 



To find the relations between the parameters of the 
family, we determined the parameters for 4 different val- 
ues of the mass Mi. The results are given in Table 4 of 
Appendix II. 

For all the members of the family, the parameters of 
Table 2 are the same and the parameters 0o and V'o are 
also the same. Only the parameters Tp, Tt and Va change 
when the mass changes. Moreover, they are correlated and 
in Table 4 (Appendix II), the relations between them are 
given. 

The precession period is proportional to the inverse of 
the square root of the mass Mi, i.e. 



Tp oc 



1/2 



(40) 



the ratio Tp/Tf, is the same for all the members of the 
family, so 



Tp « 1.967 Tb 



(41) 



and the Alfven speed is mostly proportional to the inverse 
of the precession period, i.e. 



Va OC 1/r;'-^^'* oc 



.rO.477 



(42) 



Thus, if the parameters of a member of the family are 
known, we can deduce the parameters for the rest of the 
members of the family, as Mi varies. 

From VLBI observations alone, we cannot deter- 
mine the mass of the BBH system in the nucleus of 
S5 1803-1-784. However from the knowledge of the fam- 
ily of solutions, we can illustrate the characteristics of the 
BBH system assuming a range of masses corresponding, 
for instance, to the range between the solutions Sic and 
Sid, i.e. 7.14 10^ < Mi < 5.92 10^ Mq. We find that 
the precession period is 2557.5 > Tp > 892.0 yr and the 
orbiting period is 1299.7 >Tb> 453.2 yr. 

When Ml = M2, the size of the BBH system is the 
same for all the members of the family and is 0.100 mas. 

4.7. The families of solutions with Mi ^ M2 

In the previous sections we looked for solutions assuming 
Ml = M2. 

However, it is possible to find solutions with Mi ^ M2, 
starting from the solutions found previously. To illustrate 
this possibility, we start from the solution Sic of the pre- 
vious section (see Appendix II) and find the family of so- 
lutions with Ml ^ M2 which produces the same fit. 

Assuming io — 5.8°, to — 1984.51 and rirad = 165, we 
gradually vary the mass Mi between 1.5 10^ and 1.6 10"'^^ 
and we determine M2 minimizing for each free pa- 
rameter at each step. The corresponding range of M2 is 
1.9 10^ < M2 < 1.6 10^° and the corresponding ratio 
M1/M2 is 0.01 < M1/M2 < 10. 



1e+6 1e+7 1e+8 1e+9 1e+10 1e+11 

Ml 

Fig. 9. The family of possible BBH systems with Mi 
M2 corresponding to the solution Sic of Table 4 from 
Appendix II. 



The family has been plotted in Figure [9l 
The precession period corresponding to different ratios 
M1/M2 has been plotted in Figure [TUl 




M1/M2 

Fig. 10. The precession period as a function of the ratio 
M1/M2. 

For all the members of the family solution Sic with 
Ml ^ M2, all the parameters of Table 2 are the same and 
the parameters 0o and ipo are also the same. Only the 
parameters Tp, Ti, and Va change when the ratio M1/M2 
changes. 

The changes of the parameters Tp , Tf, and Va are char- 
acterized by 

1. when Ml -> 0, M2 -> 1.98 10« Mq , 

2. when Mi < M2, Tp « Cst « 2480 yr, 

3. when Mi ^ M2, Tp increases proportionaly with 
M1/M2, 

4. Tp « 1.967 Tfa , 
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5. the Alfven speed Va is mostly proportional to the in- 



verse of Tp, i.e. 



0.958 



(43) 



6. the size of the BBH system varies from Rbbh ~ 0.050 
mas when M1/M2 — > to Rbth ~ 0.558 mas when 
M1/M2 = 10 see Figure HH 



Fig. 11. The size of the BBH system as a function of the 
ratio M1/M2. 

Therefore, the main difference between the variations 
of the parameters Tp, Tb and Va for families with Mi — M2 
and Ml ^ M2 concerns the relation between Tp and Mi. 
When Ml = M2, we have Tp oc Mj"^/^ and when Mi ^ 
M2 the variations between Tp and M1/M2 are shown in 
figure ([TO)) . 

5. Discussion and Conclusion 

Assuming that the nucleus of a radio source contains a 
BBH system, we presented a new method to fit the varia- 
tions of both coordinates of a VLBI component as a func- 
tion of time. This method is self consistent and solves the 
problems of the previous method, the two step method, 
used in Lobanov & Roland (2005). The presence of a BBH 
system produces two perturbations of the magnetic tube 
in which the plasma, that is responsible for the radio emis- 
sion of the VLBI component, propagates. These two per- 
turbations are due to the precession of the accretion disk 
and the motion of the black holes around the gravity cen- 
ter of the BBH system. The model is a geometric model 
and as we fit only the VLBI coodinates, the problem re- 
duces to an astrometric problem. As the knowledge of the 
variations of the coordinates contains kinematical infor- 
mation, we are able to deduce the inclination angle of the 
source and the bulk Lorentz factor of the ejected VLBI 
component. The fit must be done following the trajectory 
of one VLBI component and not using data of several 
components. Indeed, as a VLBI component follows the 



perturbed magnetic tube, two VLBI components will not 
follow exactly the same trajectories because the magnetic 
field perturbation propagates. In addition to the two per- 
turbations due to the BBH system, we can have a third 
one due to the slow motion of the BBH system around the 
gravity center of the galaxy. This third perturbation can 
be responsible for the slow bending of the VLBI jet that 
is often observed in compact radio sources. When it is ob- 
served, it is generally at a distance from the core larger 
than a few mas. Then studying the motion of the VLBI 
component in the first mas, we can expect that the infiu- 
ence of the slow motion of the BBH system around the 
galactic gravity center is negligible compared to the per- 
turbations due to the precession of the accretion disk and 
the motion of the black holes around the gravity center of 
the BBH system. So to obtain the characteristics of the 
BBH system in the nucleus of the radio galaxy, we have 
to fit the data concerning the VLBI component within 1 
or 2 mas only. 

To illustrate the method, we applied it to 1803-1-784 
for which we have 20 years of observations. This allows us 
to follow the VLBI component as it propagates through 
the two first perturbations of the magnetic tube. These 
perturbations correspond to CO and CI, which have been 
called oscillatory VLBI components by observers. 

Assuming, at first, that the two black holes have the 
same mass, we find for 1803+784 that 



1. the inclination of the source is io = 5.8 

2. the bulk Lorentz factor of the ejected component is 



+1.7 

-1.8 ' 



7c - 3.7 



-0.2 1 



3. the angle between the accretion disk and the plane of 
the BBH system is O = 1.09 

4. the size of the BBH system is « 0.100 mas , 

5. the precession period of the accretion disk ,Tp, and the 
orbiting period of the BBH system, T^, are related by 
Tp^ 1.967 Tfc, 

6. the origin of the component is to — 1984.51 t.^'^'l , 

7. the duration of the ejection of the plasma responsible 
for the VLBI component is {Tejec)obs = 0.93 to'34 

the observer frame which corresponds to {Tejec)bbh ~ 
13.9 yr in the BBH frame. 



We find that there is not a unique solution for the 
mass of the BBH system, but a family of solutions that 
produces the same fit. For all the members of the family, 
all the parameters are the same, except the parameters 
Tp, Tfc and Va which change when the mass changes. The 
variation laws for these parameters are : 



Tp oc 



-1/2 



Va OC l/Tp"-954 oc Mi"-^^^ 



(44) 



(45) 



the ratio Tp/Tt is the same for all the members of the 
family, so 



Tp « 1.967 Tt , 



(46) 
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Thus, if the parameters of a member of the familly are 
calculated, we can subsequently deduce the parameters of 
all the members of the family when Mi varies. 

To illustrate the characteristics of the BBH system, if 
we assume a range of masses corresponding to 7.14 10^ 
Mq < Ml < 5.92 10^ Mq, we find that the precession 
period is 2557.5 yr > Tp > 892.0 yr, the orbit period is 
1299.7 yr >Tb> 453.2 yr . 

Furthermore, from the knowledge of the solution with 
Ml = M2, we can find the families with Mi ^ M2- 

If VLBI observations allow us to find the possible fam- 
ilies of BBH system solutions and the main characteristics 
of these systems, radio sources exist for which, in addition 
to VLBI data, optical bursts associated with the birth of 
the VLBI component have been observed (see for instance 
Britzen et al. 2001 in the case of 0420-014 and Lobanov & 
Roland 2005 for 3C 345). The combination of VLBI and 
optical observations can further constrain the families of 
solutions. 
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6. Appendix I : CoefFicients A, B and C 

Let us call 



4){t) = ujpt - kpz{t) + 4>o 



and 



tjj{t) = LJbt - kbz{t) + t/^o 



With 



Mo 



xi = yi 



-,1/3 



(47) 



(48) 



(49) 



Ml + M2 

the coefBcients A, B and C of equation ([^ are given by 



(xi cos(V'o))^ + {yi sin(V'o))^ 



7^2 



C = eM-^t/Td) 

/ yicjj,yi sin('0o) ^ xia;i cos(V'o) 
I Td Tj 
f xiuJbXiCoa{ipo) yiyism{i>o) 



2R{z) 



2R{z) 



yi sin(V'o)wp xi cos^ipo) 



Td 



^d 



xi cos{ipo)ujp yi sm{^o) 



Td 



^d 



cos{tp{t)) 
sm{^{t)) - 
cos((/)(t)) 
sin(0(i)) - 



sm{ij{t) - cb{t)) {-{iJb - ojp)R{z) {xi + yi) /Td} + 
cos('0(t) - (t>{t)) {ujpLUbRiz) {xi + yi) - 

R{z) (x, + y,) /Tj] + '^{xl+ yl) + R\z)ul + 



R^{z) + 



xl+yl 



/Tj 



(52) 



A = exp(-2i/Td) 



LOb^pR{'- 



-{yi+xi)cos{xJj{t)-^{t)y 



^^(a^i + yi) sin(V(t) - <^(i)) + 



(50) 



exp(-2t/rd) 
2yiUJbyi sin(V'o 



2a;iW6a;i cos(i/'o) . , , 

-sin(i/'(0)- 



di?(z) 



dz 
f dR{z) 
\ dz 



Xi C0s(V'o 



yi sin(V'o) 



cos(7/;(t)) + 

yi sm{'4)o)Riz)u;p 



Va 

xi cos{ipo)R{z)ujp 

Va 



cos{(j){t)) 
sin(0(t)) - 



sin(V'(i) - 0(t)) <i -^b^^ {xi + yi) 



{xi+yi)/Td} + 



{ujp - ujb)R{z) 

Va 

COs(^(t)-0(t))(- '"^"^^^") (XI +,1)- 
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dz ^ — "J 
2ljIR^{z) 2dR{z)R{z) 



+ yl) 



Va 



dz Td 



(51) 



Roland & al.: Modelisation of compact radio 

7. Appendix II : BBH parameters for 4 values of 
Ml of the family SI 

The BBH parameters for 4 values of Mi of the family SI 
are given in Table 4. 



Table 4 : BBH parameters for 4 values of Mi of the 

family SI 





Sla 


Sib 


Sic 


5*1^ 




78.9 


78.8 


79.0 


79.1 




16380.4 


6608.55 


2557.47 


891.95 


Ml 


1.7556 10^ 


1.0743 10** 


7.1437 10** 


5.9232 10^ 




59.5 


59.5 


59.5 


59.8 


n 


8326.9 


3360.29 


1299.68 


453.198 


Va 


0.0078855 


0.019298 


0.048358 


0.12665 



For all the m,em,bers of the family, all the parameters 
of Table 1 are the same and the parameters (po and ■i/'o are 
also the same. Only the parameters Tp, and Va change 
when the mass Mi changes. 



